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We reexamine the numerical study of the magnetized state 
of the S = 3/2 spin chain with single ion anisotropy D(> 0) 
for the magnetization M = Ms/3, where Ms is the satura- 
tion magnetization. We find at this magnetization that for 
D < D c i — 0.387 the system is critical and the magnetiza- 
tion plateau does not appear. For D > D c i, the parameter 
region is divided into two parts D c \ < D < D C 2 = 0.943 
and Dc2 < D. In each region, the system is gapful and the 
M — Ms/3 magnetization plateau appears in the magneti- 
zation process. From our numerical calculation, the interme- 
diate region D c i < D < D C 2 should be characterized by a 
magnetized valence-bond-solid state. 
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I. INTRODUCTION 

Recently there has been a considerable interest in the 
magnetization process of one-dimensional quantum spin, 
systems. Extending the Lieb-Schultz-Mattis theorem,El 
Oshikawa et aln gave a necessary condition for the ap- 
pearance of the magnetization plateau as p(S — (m)) = 
integer, where S is the magnitude of the spin, p is the 
periodicity of the magnetic ground state in the thermo- 
dynamic limit, and (m) is the magnetization per site. 

As a simple case, the 5 = 3/2 spin chain with single 
ion anisotropy in a magnetic filed h, 
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has been studied. Here Sj is the 5 = 3/2 spin at the j-th 
site, and M is the magnetization M — J^j Sj- Hereafter, 
we define the saturation magnetization as Ms(= 3L/2). 
In this model, a magnetization plateau can appear at 
M = M s /3 (p = l). For D = (S = 3/2 Heiseru 
berg model), a magnetization plateau does not appear n 
On the other hand for large D, the spin of each site 
tend to have S z = 1/2 for M — Ms/3, and a mag- 
netization plateau appears. Oshikawa et alia checked 
that a magnetisation plateau appears for D = 2. Sakai 
and Takahashia studied the model with the numerical 
diagonalization calculation using the phenomenological 
renormalization group (PRG) analysis, and found that 
the transition point is D c = 0.93 ± 0.01. They con- 
cluded this point_aa. the Berezinskii-Kosterlitz-Thouless 
(BKT) transitionEm point, that is, the plateau and the 



no-plateau transition point. Although the PRG analysis 
is a powerful method to study the second order phase 
transition, this is not sufficient for the BKT transition 
because of the singular behavior of the energy gap and 
of the logarithmic size dependence of excitation energies. 
Thus we think that it is worthwhile to reexamine this 
problem by use of another numerical approach. 

The organization of this paper is as follows. In the 
next section, we consider an effective continuum model 
in order to investigate the mechanism of the appearance 
of the M = M5/3 magnetization. In section 3, we show 
the numerical method and the results with the numeri- 
cal diagonalization data. To study the system, we firstly 
consider the finite-size behavior with conformal field the- 
ory and the renormalization group analysis. Next, using 
them, we determine the boundary of the plateau and the 
no-plateau regions, and also check the critical behavior. 
Last section is devoted to summary and discussion. 



II. EFFECTIVE CONTINUUM MODEL 

For the magnetization M = Ms/3, the effective con- 
tinuum model describing the low energy physics is given 
byH 
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where n is the momentum density conjugate to 0, 
[c/>(x), n(aj')] = iS(x — x'), vs is the sound velocity of 
the system. The dual field 9 is defined as d x 9 — nil. The 
correlation exponents of several operators are governed 
by the coupling K and the scaling dimension of the op- 
erator v^cosv^f^ is x — K/2. Thus the second term in 
the Hamiltonian (0) is relevant for K < 4 and irrelevant 
for K > 4. 

The one-loop renormalization group equations of this 
model are given by 



dK(L)- 
dlnL 
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where L is an infrared, cutoff. These are the recursion 
relations of Kosterlitz.13 The renormalization group flow 
is shown in Figure 1. In the region A, the coupling K 
renormalized to a finite value and y^ to zero; the system 
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is critical and this corresponds to the no-plateau region in 
the spin system. In the regions B and B', K renormalized 
to and \y^\ to oo; these regions are massive phases and 
the field <fi is locked to tt and for B and B' respectively. 
There is a separatrix 32K~ 1 - 8 In if -1 - = 8 + 81n4 
which separates the infrared stable region (A) from the 
infrared unstable regions (B and B'). On this separatrix 
the BKT transition takes place. Near (K,y,p) = (4,0), 
and with the notation K = 4(1 + yo/2), the BKT transi- 
tion point is on the line jy^j — yo and we obtain 
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where yo and L$ are the bare values. 

According to the numerical calculation by Oshikawa 
et alB for the S = 3/2 Heisenberg chain (D = 0) at the 
magnetization M — Ms/3, this Gaussian coupling K is 
K 4.4, (in the notation of ref.0 the compactification 
radius is given by R = ^ttK)- 1 / 2 = 0.95/^8^), and 
the system is gapless (no plateau). This value is slightly 
larger than the BKT critical (fixed point) coupling K = 
4. 



III. NUMERICAL APPROACH AND RESULTS 

A. Finite size behavior 

In order to study numerically, let us consider the finite- 
size behavior of the effective model (||) based on the con- 
formal invarianceJa The scaling dimension of the primary 
field 0„ hn — exp(im\/20 + in\/29) for the fixed point 
Hamiltonian (y^ — 0) with the periodic boundary condi- 
tion (PBC) is given by 
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where m and n are the "magnetic" and the "electric" 
charges in the Coulomb gas.EI According to the confor- 
mal field theory, the state of the finite-size system have 
a correspondence to the operator in the infinite-size sys- 
tem, and the excitation energy of the finite-size system 
is given by 



AE m ^ n (L) — E„ l:n (L) — E g (L) = 
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where E g (L) is the ground state energy of L-site with 
PBC. 

The BKT transition takes place between the gapless 
and gapful regions, and the PRG analysis is affected by 
several irrelevant corrections sensitively. Moreover, on 
the transition point, there exists a logarithmic size cor- 
rection in the excitation energy. PRG analysis does not 
suppose these features. In order to determine the BKT 
transition point, we .apply the methpd proposed by No- 
mura and Kitazawat3 (see also retail andEJ), which is 



free from the logarithmic size correction. If we can have 
the half-integer for the "magnetic" charge m, which is 
not permitted for PBC by the translational symmetry, 
we have the following finite size spectrum for m = ±1/2, 
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where the energies AE^y 2 Q (L) correspond to the op- 
erators -\/2cos</>/-\/2 and V2 sin</>/\/2 respectively. The 
term with y^ comes from the first order-perturbation and 
from the operator product expansionsjLj 
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Near the BKT transition point \y^,\ = yo (K = 4(1 + 
yo/2)), defining = ±yo(l + t) where t measures the 
distance from the transition point, we have for > 



c ±i/2,o(-k) — 2 
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and for < 0, 
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On the other hand, for (m, n) = (0, ±2) we have 
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Xo < ±2(L) = 2^r s AEo ' ±2{L) = i = \- i yo{L) - (9) 

From cqs. (Q), (g) and (g), energy differences 
AE±y 2 (L) and AEo t ±i cross linearly at the BKT tran- 
sition point t = for y<p > (y<j, < 0), and this behavior 
can be used to determine the BKT critical point. ■— ■ 
From (||), AE^_ 1 ^ 20 and AE s ±l ^ 20 cross at = 0Hl 
This is the Gaussian fixed point. For K < 4, the operator 
V2 cos \/2<f>, which is in the second term of the Hamilto- 
nian (Q), is relevant and a second order phase transition 
occurs at this point. On this transition point, we have 

AE ±l/2,0 < AE 0,±2- 

The "electric charge" n in eq. (||) relates to the vari- 
ation of the magnetization from M — Mg/3, and the 
excitation energy Ai?o,±2 in the spin system is described 

as 



A£ ,± 2 (£) 
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where E(M, L) is the lowest energy for the magnetization 
M with the PBC. The excitation energies corresponding 
to the operator v^cos^/v^ and \/2 sin <p / \[2 are ob- 
tained by the two lowest-energies with the twisted bound- 
ary condition (TBC)ll3~E] sf = -Sf, S Z L+1 = Sf 



as 



AE c ±1/2fi (L) = E l ^(M s /3,L,l) - E(M S /3,L) (11) 
A^iiAoW = E TBC (M S /3,L,-1) -E(M S /3,L) 

where E TBC '(M s /3, L, P) is the lowest energy with the 
magnetization M = Ms/3 and the parity P : Sj — > 



B. Numerical results 

Here we show our numerical results. We calculate the 
energy of finite-size systems L = 8, 10, 12, 14 with the nu- 
merical diagonalization calculation. In the region D > 0, 
the lowest energy state with the magnetization M and 
with PBC has the parity P = (-1) M + L / 2 and the mo- 
mentum q = (M + L/2) x 7T. 

Figure 2 shows the energy differences A£|j, 2 Q , 
A£|j, 20 , and AP ,±2 for L = 14 systems. We find 
that among these three energies the lowest one is APo,±2 
for small D (0 < D < D Cl ) and A£^ 1/2 Q for large Z? 
(D C 2 < D). In the intermediate region D c \ < D < D C 2, 
Ai?J 1 y 2 is lower than AP^. 1 ^ 2 and AP ,±2- In Figure 
1, the small D region < D < D c \ should correspond to 
the region A, the intermediate region D c \ < D < D C 2 to 
B, and the large D region to B\ 

From the critical properties of the effective model (Q) 
[see eqs. (||) and (9)], the point D c \ should be a BKT 
critical point. The size dependence of the crossing points 
between A£J. 1(/2 and A£'o.±2 is shown in Figure 3 and 
the extrapolated value is given by D c \ — 0.387. For small 
D (< D c i) the excitation spectrum is gapless and there 
does not exist a plateau in the magnetization curve. 

In order to check this, we calculate the averaged scaling 
dimension,!^ E3 
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at the point D = D c \. From (0), this value cancels the 
leading logarithmic size correction, and should be 1/2 at 
D = D c i. To calculate the scaling dimension, we need 
the sound velocity ug. This can be calculated using the 
lowest energy with M — Ms/3 and the wave-number 
q = 2-k/L (corresponding to the U(l) current), 
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We extrapolate this finite size value in the polynomial 
of l/L 2 , and obtain 2-kvs — 20.06. Figure 4 shows the 
size dependence of the value (|l2]). We can see that the 



extrapolated value is very close to the expected one 1/2. 
We also calculate the conformal anomaly number c from 
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where e m=1 /2 is the energy per site of the limit L — > oo, 
c is the conformal anomaly number, and • • • means the 
higher size correction. We estimate c at D = D c i as 
c = 0.97, which is close to the expected value c = 1. 

Thus we can conclude that the point D = _D cl is 
the BKT critical point. The width of the plateau 
near the transition point behaves singularly as Ah ~ 
exp ( — const . / y/D — D c \) . 

According to (g), the crossing point D C 2 between 
AP^_ 1 ^ 2 Q and APJ 1 ^ 2 should be a second order phase 
transition point (K < 4 and = 0). For large D 
(D C 2 < D), the system is in the large- D region and 
a plateau appears in the magnetization curve at M = 
Ms I '3. Thus we have another massive phase in the re- 
gion D c i < D < D C 2- As a possibility, this intermediate 
region should be characterized by the following partially 
magnetized valence-bond-solid (VBS) state, 
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where we describe the spin state by the Schwinger bosons, 
that is, a] (&t) creates the S = 1/2 t (1) spin at the j- 
th site, and |0) is the vacuum of the Schwinger bosons. 
Here we assume PBC. This partially magnetised VBS 
state had firstly been discussed by OshikawaEll for an- 
other model. (see also ref.13) 

Let us consider the crossing behavior near the point 
D = D C 2 on another point of view. The large- D phase is 
characterized by the state 



|large-D> = nat(aj&t)|0>. 



(16) 
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With the twisted boundary condition (a L , t — a\, 

—b\), the magnetized valence-bond-solid state is de- 
scribed as 
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while the large-D state is also characterized by ( |16| ) . The 
state |mVBS,TBC) has the parity of the space inversion 
P = — 1, while the large- D state has P = 1. This ex- 
plains the level crossing between AB^wj Q and AE^^ Q 
in Figure 2. 

The size dependence of the crossing point is shown in 
Figure 5. The estimated critical point between the inter- 
mediate and the large- D regions is D c = 0.943, which is 
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very close to the value D c = 0.93±0.01 obtained by Sakai 
and Takahashi. However, the phase transition is not of 
the BKT type but a second order phase transition which 
is the same universality class as the Haldanc,- large-D 
transition in the 5=1 anisotropic model, E-J1I3 

L L 

H 3=1 = £[5J5* +1 + S]S* +1 + A5|5| +1 ] + DY+s*f. 

3=1 3 = 1 

(18) 

Especially, the intermediate region of the model ([!]) cor- 
responds to the 5=1 Haldane phase in (|l^). 

To verify this speculation, we calculate the coupling K. 
The sound velocity is evaluated as 2ttvs = 21.87. Using 
this value, we determine K from x^_ 1 y 2 (L), x±^, 2 (I) 
[eqs. (0)1, and 
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as shown in Figure 6. Extrapolated values of these three 
are consistent and we obtain K = 3.77. This is smaller 
than the BKT fixed point value, and the point D C 2 = 
0.943 should be a second order phase transition point. 
We also estimate the conformal anomaly number form 
cq. (14), and check that it is very close to the expected 
value c = 1. 

In the vicinity of the critical point D — D C 2, the width 
of the magnetization plateau decreases as Ah ~ \D — 
D C 2\ l where the critical exponent v is given by v — 
1/ (2 - x± h0 ) = 8.68 



mapped the 5 = 3/2 model to an S = 1 generalized 
anisotropic model with h = 0. In order to investigate the 
relation between the 5 = 1 with h = and the 5 = 3/2 
model in a magnetic field, we further applied the same 
analysis to the following model, 

#A = J2(SjS'Ul + S 3 S j+l + AS 3 S 3 + ^ 

3 

L L 

+ D Y,( s D 2 - h Y, s h (2°) 

3 3 

in the region < A < 1. Figure 7 shows the energy 
differences AE± 1/2 Q , AE s ±1/2 , AE ,±2 for L = 14, A = 
0.5 system. The lowest energy difference is Ai?o,±2 for 
small D and AE^,^ for large D. In this case, there 
is a transition between the no-plateau and the large- 1? 
regions. [The crossing point of AE^_ ± ^ 2 Q and AE^ , 2 
is the Gaussian fixed point = in the no-plateau 
region (K > 4).] 

We show the A-D diagram in Figure 8. We found that 
for 0.808 < A < 1, the M = M s /3 phase structure is 
the same as is of the A = 1 case and there exists the 
intermediate magnetic region (D c i < D < D^). But for 
< A < 0.808, the direct transition (BKT type) occurs 
between the no-plateau and the large-D regions. This 
boundary runs from (A, D) = (0, 1.312) to the multicrit- 
ical point (0.808, 0.817). The phase structure is topplop^ 
ically the same as is of the 5=1 anisotropic modelE-3~tL£l 
with h = (U). 

We acknowledge K. Nomura, M. Oshikawa, and T. 
Sakai for discussions. The computation in this work was 
performed using the facilities of the Supercomputer Cen- 
ter, Institute for Solid State Physics, University of Tokyo. 



IV. SUMMARY AND DISCUSSION 

We studied the 5 = 3/2 spin chain with the single ion 
anisotropy in a magnetic field (|l]). With the numerical 
approach based on the field theoretical analysis, we found 
a new region between the no-plateau (small D) and the 
magnetized large- 1? (large D) regions. As far as we know, 
this is the first work in which no-plateau state and two 
kinds of plateau states are found when the parameter is 
changed in a simple and realistic model. This intermedi- 
ate region should have the same character of the partially 
magnetized valence-bond-solid state, and-sbeiuld be char- 
acterized by the string order parametercilia-and by the 
edge states for the open boupdacy systemEJcJ reflecting 
a hidden discrete symmetryHrEa But in our small size 
calculation, we could not collect evidence, which remains 
for future study. 

Restricting to three states S z — 3/2, 1/2, and —1/2, 
(in sufficiently large magnetic field) Sakai and Takahashi 
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FIG. 1. Renormalization group flow diagram. The bold 
lines (separating the regions A and B, and A and B') are the 
BKT critical lines. On the dotted line between the regions 
B and B' (y^ = 0, K < 4), a second order phase transition 
occurs. 
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FIG. 2. Energy differences A£Jy 2i0 {o),AE s ±1/2>{) (.),and 
A.Eo,±i (O) for the L = 14 system. The crossing point of 
AByyj and ASo,±2 is the BKT critical point. The cross- 
ing point of and A£|_ 1 / 2 is a second order phase 
transition point. 
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FIG. 3. Size dependence of the crossing point between the 
energy differences AE± 1 / 2 and A_Eo,±2- The BKT critical 
point is estimated as D c \ = 0.387. 
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FIG. 4. Size dependence of the averaged scaling dimension 
(^±1/2,0 + 3:r ±i/2,o)/ 4 at the BKT critical point D cl = 0.387. 
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FIG. 5. Size dependence of the crossing point between the 
energy differences AE± 1 / 2 and A£|. 1 y 2 . The critical point 
between the intermediate and the magnetic large-D regions 
is estimated as D C 2 = 0.943. 
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FIG. 8. The boundary of the no-plateau, the large-D, and 
the magnetized VBS (m-v) regions. The boundary between 
the no-plateau and the large-D regions, and between the 
no-plateau and the magnetized VBS regions (o) are the BKT 
critical lines. The boundary between the large-D and the 
magnetized VBS regions (O) is a second order critical line. 
The multicritical point of these regions (•) are estimated as 
(A, D) = (0.808, 0.817). We also show the Gaussian fixed line 
in the no- plateau region (x). 
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FIG. 6. The parameter K obtained by the data 

*±l/2,0 = K / 8 (°)> *±l/2,0 = K / 8 W. aIld ^,±1 = 1/2K 

(O). The extrapolated value is K = 3.77. 
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FIG. 7. Energy differences A£^ 1/2 (o), AE± 1/2 (•), and 
A-E ,±2 (O) for the system with A = 0.5 and L = 14. The 
crossing point of AEj. 1 / 2 ,o and AEo,±2 is the boundary be- 
tween the no-plateau and the large-D regions. The crossing 
point of AE^^ g and AE s ±1 / 2 is the Gaussian fixed point 
in the no-plateau region. 
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